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Abstract 

We observed spectroscopically the eclipsing binary system T-Cygl-01385 in 
order to determine physical properties of the components. The double-lined 
nature of the system is revealed for the first time and the radial velocities 
are obtained for both stars. We have derived masses, radii and luminosi- 
ties for both components. Analyses of the radial velocities and the Kepler- 
Cam and the TrES light curves yielded masses of Mi = 1.059±0.032 M & and 
M 2 =0.342±0.017 M and radii of R!=1.989±0.022 R Q and R 2 =0.457±0.013 
Rq. Locations of the low-mass companion in the mass-radius and mass- 
effective temperature planes and comparison with the other low-mass stars 
show that the secondary star appears just at the transition from partially to 
fully convective interiors for the M dwarfs. When compared to stellar evolu- 
tion models, the luminosities and effective temperatures of the components 
are consistent with Z=0.004 and an age of about 6 Gyr. A distance to the 
system was calculated as d=355±7 pc using the BV and JHK magnitudes. 
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method: spectroscopy 



1. Introduction 



The lower-mass stars of the main-sequence of the Hertzsprung-Russell 
diagram are known as M-dwarfs. They constitute the majority of stars in 
the solar neighbourhood. They are intrinsically faint since they are cooler 
and smaller than the other stars . Because of their faintness photometric and 
spectrosc o pic o 



Dservations could be obtained for a limited number M-dwarfs. 
Chabrierl (120031 ) suggests that at least 70 per cent of all stars in the sphere 
with a radius of 10 pc about the Sun are M-dwarfs. Therefore, detailed 
studies of these faint but numerous low-mass stars are critical importance to 
evaluate a model of the galactic evolution and present status. Moreover the 
M-dwarfs are particularly needed in the understanding of evolution of main- 
sequence stars towards the lower-mass regime. The evolutionary calculations 
of low-mass stars allow us to define not only the transition from partly to 
fully convection b ut also to define a lim it between stars and brown dwarfs. 



As pointed out by 



Shulyak et al. 



( 1201 ll ) around spectral type M3.5 stars be- 
come fully convective and thus the dynamo mechanism must be different in 
cooler stars because they do not posses ta c hoclin e layer with strong differ- 



ential rotation. However, 



Reiners fc Basril ( 120071 ) have concluded from the 



measurements of Stokes parameter about tens of stars spanning the whole M- 
dwarfs that no significant change in the average magnetic field strengths oc- 



curs w hen stars become fully convective. Later on lReiners. Basri.fc Browning 
( 120091 ) could reveal the rotation- magnetic field relation in which magnetic 
field strengths increase towards short rotation periods. 
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Low-mass M-dwarfs are generally faint objects because of their small radii 
and low temperatures. Due to their faintness the M-dwarfs have limited the 
number of high-resolution, high signal-to-noise spectroscopic studies. In ad- 
dition, the optical spectrum of these cool stars is mostly covered by molec- 
ular bands which hide and blend the atomic lines used in spectral analysis. 
These dominant molecular bands make it difficult to measure the atomic line 
strengths which are needed for metallicity determination. The physical pa- 
rameters, such as mass, radius and luminosity as well as age and rotational 
period could only be derived if a M-dwarf is a member of a close binary or a 
multiple system. In a binary or a multiple system the components are most 
likely coeval and their spin axes are perpendicular to the orbital plane. 

Numerous photometric observat ions of binary stars were re cently gath- 



ered by the surveys of like the TrES 



2004 ). SuperWASP flChristian et al. 



Alonso et al. 



1996 ). NSVS (IWozniak et al. 



20061 ) , Kepler (IBorucki. Koch. &: Kepler Science Team! . 



20101 ) . etc. Despite the main aim of these surveys is to search gamma- 



ray bursters, and especially extra-solar planetary transits many binary sys- 
tems are discovered and light curves of many systems could be obtained. 
Therefore, the tremendous photometric datasets containing unknown or lit- 
tle known binary systems are presented for the use of ast ronome rs. 



Light va r iabilit y of T-Cygl-01385 was announced by iDevorl (120081 ) and 



Devor et al. 



(120081 ) in the list of 773 eclipsing binaries found in the Trans- 
Atlantic Exoplanet Survey. It was classified as an "ambiguous" binary in 
their lis t. Mass and rad i us of the components were estimated for the first 



time by 



Fernandez et al. 



(120091 ) combining the spectroscopic orbital elements 



obtained from the primary star's radial velocities with a high-precision transit 
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Figure 1: The residuals between the observed times of minimum light and computed with 
the new ephcmcris. 

light curve obtained by the KeplerCam. The results obtained up-to-date 
point out that the fainter component of T-Cygl-01385 is so close to the fully 
convective M-dwarfs. Therefore we planned new spectroscopic observations 
of the system to refine the masses, radii and effective temperatures of both 
stars. 



2. Period determination 

The catalog information for T-Cygl-01385 was given in Table 1. First 
we collected the times for mid-light minimum obtained by various automatic 
and robotized telescopes and surveys. These times of minima are given in 
Table 2 as averaged for the filters used. A linear least-squares fit to the data 
yields the following ephemeris 

Mini = HJD 2453926.7896(3) + 6 d . 558662(31) x E 

where the bracketed quantity is the uncertainty in the last digit (s) of the 
preceding number. All uncertainties quoted in this work are standard errors. 
The residuals of the fit are plotted in Fig. 1 and show no indication of any 
form of period change in about ten years. 
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Table 1: Catalog information for T-Cygl-01385. 
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2008), C GSC2.3 G uide S tar Catalog, version 2.3.2 (Morrison ct al 



Perrvman fc ESA , 



1997), e NOMAD NOMAD Catal og, l|Zacharias et al 



1999) 



20051 ). 1 CMC14 Carlsberg Meridian Catalog 14 Jcopenhagen University et allbooa ). 9 ASCC-2.5V3 All- 
Sky Compiled Catalog of 2.5 Millon Stars, jKharchenko fc Roesel I2OO9I) , h ACT A strographic Catalog , 



I Urban. Corbin, fc Wvcoa , ll997l) . i TASS The Amateur Sky S urvey (TASS) Catal og, foroege et al 



2006), 



i UCAC High den sity, Highly Accura te, Astrometric Catalog, llUrban et al.Ll2004T) . * 7V E g Trans Eclipsi n. 



Table 2: Times of minimum light for T-Cygl-01385. The 0-C values refer to the difference 
between the observed and calculated times of mid-eclipse. 
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Ref: 



fi) — ([Copenh age n Uii iversi 



Fernandez et al 



mm, (4) 



v et al 



Christian et al 



(|2006|) 



200C)j (S) iDevori (|2008h . (3) 



3. Echelle Spectroscopy 

Optical spectroscopic observations of the T-Cygl-01385 were obtained 
with the RTT15C0 1.5 m telescope at TUg[] using the R ~ 7000 Echelle 
spectrograph that covers 400 nm < A < 900 nm. These observations were 
used to resolve the components of T-Cygl-01385 and get individual radial 
velocity measurements for each star in the system. Eleven spectra were 
obtained using 60 minute integrations on 8 nights in July, 2010 with typical 
signal-to-noise ratios of ~120 at 6563 A. 

The spectra were processed in the standard way for cross-dispersed Echelle 

1 http://tug. tug. tubitak.gov.tr/rttl50_tfosc.php 
2 TUBITAK National Observatory (Turkey) 
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spectra, using FXCOR package in IRAF. The routine processes the data using 
biases and halogen lamp observations taken at the beginning of the night, 
median combines three individual images while performing cosmic ray rejec- 
tion, extracts the individual orders from the combined image and performs 
the wavelength solution on each order using a FeArarc lamp taken either 
before or after each set of science exposures. 

A total of 11 orders are used each night to derive radial velocities via 
cross-correlation with a standard template. We use the bright tPsc (F7V) 
and 50 Ser (F0 V) as the heliocentric radial velocity standard stars. Each 
spectral order is cross-correlated separately, then an iterative 3-er clipping 
is performed prior to performing a weighted average to obtain a final radial 
velocity measurement for each night. The components are identified each 
night via the peak and width of each feature in the cross correlated functions. 
The typical radial velocity precisions ranged from 2 to 13 km s -1 for the 
various components. 

3.1. Radial velocities 

To derive the radial velocities, the eleven spectra obtained for the system 
are cross-correlated against the template spectra of st andard stars on an 



order- by-order basis using the FXCOR package in IRAF ( jSimkinl . I1974J ). 

The majority of the spectra showed two distinct cross-correlation peaks 
in the quadratures, one for each component of the binary. Thus, both peaks 
are fitted independently with a Gaussian profile to measure the velocities 
and their errors for the individual components. If the two peaks appear 
blended, a double Gaussian was applied to the combined profile using de- 
blend function in the task. For each of the eleven observations we then 
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determined a weighted-average radial velocity for each star from all orders 
without significant contamination by telluric absorption features. Here we 
used as weights the inverse of the variance of the radial velocity measurements 
in each order, as reported by FXCOR. 

The heliocentric radial velocities for the primary (V p ) and the secondary 
(V s ) components are listed in Table 3 , along with the dates of observa- 
tions and the corresponding orbital phases computed with the new ephemeris 
given in the previous section. The velocities in this table have been cor- 
rected to the heliocentric reference system by adopting a radial velocity value 
for the template stars. The radial velocities are plotted against the orbital 
phase in Fig. 2. The ra dial velocities of the more massive star measured by 



Fernandez et al. 



( 120091 ) are also plotted as circles in the same figure. There 
is no systematic difference between the measurements. 

We analysed the radial velociti es for the initial orbita l para meters us- 



ing the RVSIM software program (IKane. Schneider, fc Gel . 120071 ) . Figure 2 
shows the best-fit orbital solution to the radial velocity data. The results 
of the analysis follow: 7= -9±1 km s" 1 , #i=31±3 and K 2 =96±A 

km s _1 with circular orbit. Using these values we estimate the projected or- 
bital semi-major axis and mass ratio as: asim=16.46±0.18 R@ and q = 
f|=0.323±0.020. 



4. Light curves and their analyses 



Photometric observations of T-cygl-01385 were obtained by automatic 



and robotised telescopes. The first complete 



the TrES wide-angle transiting planet survey (lAlonso et al- 



ight curve was obta ined by 



19961 ). Addi- 
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Figure 2: Radial velocities for the components folded on a period of 6.5587 days. Symbols 
with error bars (error bars are masked by the symbol size in some cases) show the radial 
velocity measurements for the components o f the system (primary: fi lled circles, secondary: 



open squares). The velocities measured by 
circles. 



Fernandez et al 



([20091 ) are shown by empty 
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Table 3: Heliocentric radial velocities of T-Cygl-01385. 
The columns give the heliocentric Julian date, the or- 
bital phase (according to the ephemeris in §2), the radial 
velocities of the two components with the corresponding 
standard deviations. 
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Figure 3: Unbinned phase-folded light curve of T-Cygl-01385 obtained by the TrES in the 
Sloan r-passband and the best-fit model. The very shallow secondary eclipse is detectable 
at around phase 0.5. 



tional phot ometric data, especia 



the NSVS flWozniak et al. 



ly in the prima ry eclipse, were obtain ed by 



200 4|), SuperWASP 



Christian et al 



ifep/erCam (jBorucki. Koch, fc Kepler Science Team 



2006|), and 



20101 ). Since the data 



gathered by the SuperWASP have too large scatters we do not include them 
into the analysis for the orbital parameters. In Fig. 3 we plotted the TrES 
data against orbital phase. The observations obtained by the TrES and the 
KeplerCam within the primary eclipse are shown in Fig.4. 

We may constrain the effective temperature and spectral type of the pri- 
mary star using the By, Vt, J, H, and K magnitudes which are already given 
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in Table 1. We have derived V= 10.96, B-V=0.56, J-H=0.241 and H-K=0.08 



m ag. Comparing the co 



by 



Drilling fc Landolt 



or in dices with color-spectra l type calibrations given 



f l2000h and 



Tokunaga 



(120001 ) we estimated a spectral 



type of F8 sub-giant for the primary star. Thus, an effective temperature of 
T e ff = 6 250 ± 100K and a color excess of 0.04 mag are estimated. 

We used the most recent version of the ec lipsing binary light curve model- 



ing algorithm of 



Wilson fc Devinney 



(119711 ) , as implemented in the PHOEBE 



code of Prsa & Zwitter (2005). In order to obtain the physical parameters of 
the component stars we, first analysed the TrES data. The code needs some 
input parameters, which depend upon the physical structures of the com- 
ponent stars. The values of these parameters can be estimated from global 



stellar prop erties. Ther e 
cients from van Hamme 



'ore, w e adopted the linear limb-dar 



rening coeffi- 



19931 ). the bolometric albedos from 



Lucy 



(1967) 



and the gravity brightening coefficients as 0.32 for both components. The 
rotational velocities of the components are assumed to be synchronous with 
the orbital one. 

The adjustable parameters in the light curves fitting were the orbital in- 
clination, the surface potentials of the two stars, the effective temperature 
of the secondary, and the color- dependent luminosity of the hotter star, the 
zero-epoch offset, semi-major axis of the orbit, the mass-ratio and the sys- 
temic velocity. A detached configuration (Mode 2) with coupling between 
luminosity and temperature was used for solution. The iterations were car- 
ried out automatically until convergence, and a solution was defined as the 
set of parameters for which the differential corrections were smaller than the 
probable errors. Our final results are listed in Table 4. The uncertainties 
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Figure 4: Eclipse light curves of T-Cygl-01385 obtained by the TrES at the Sloanr- 
passband and the KeplerC&m at the z-passband. The continuous lines show the best 
fits. 



13 



Table 4: Results of the simultaneous analyses of the TrES and the KeplerCam light curves 
for T-Cygl-01385. 

Parameters Adopted 



i° 86.4±0.02 

T e//l (K) 6 250[Fix] 

Teffa (K) 2 940±120 

Qi 8.609±0.023 

Q 2 12.989±0.030 

ri 0.1206±0.0003 

r 2 0.0272±0.0001 

, T L 2 T ; (TrES-r) 0.9993±0.0005 

(Ll L 2 L2) (Kepler-z) 0.9996±0.0004 

X 2 0.051 



assigned to the adjusted parameters are the internal errors provided directly 
by the code. The computed light curve corresponding to the simultaneous 
light- velocity solution is compared with the observations in the Fig. 3 and 4. 

The fundamental stellar parameters for the components such as masses, 
radii, luminosities were calculated and listed in Table 5 together with their 
formal standard deviations. The standard deviations of the parameters have 
been determined by JKTABSDFvl] code, which calculates distance and other 
physical parameters usin g several different sources of bolometric corrections 



( jSouthworth et all 120051 ). The mass for the primary of Ma = 1-06 ± 0.03M 



3 and secondary of Mb = 0.34 ± 0.02M^ ar e cons isting of an evolved late 







F-star and mid M-dwarf ([Drilling fc Landolt 1 . 12000| ). 



There is no measured trigonometric parallax available for the system. 
From the B- and V-passband measurements of Tycho and the JHK magni- 



3 This can be obtained from http://http://www.astro. keele.ac.uk/~jkt/codes. html 
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Parameter 



Table 5: Fundamental parameters of T-Cygl-01385 

Primary Secondary 



Spectral Type 

a (R ) 

V 1 (km s" 1 ) 

«(°) 

g 

Mass (M Q ) 
Radius (R©) 

Teff (K) 

log (L/Lq) 
log 9 {cgs) 

(vsin i) ca i c , (km s _1 ) 
d (pc) 



F8(±2)IV-V M4(±1)V 
16.49±0.18 

-9±1 
86.36±0.02 
0.323±0.020 
1.059±0.Q32 0.342±0.017 
1.989±0.022 0.457±0.013 
6 250±100 2 940±120 
0.736±0.035 2.039±0.083 
3.866±0.005 4.652±0.026 
15±1 4±1 
355±7 



tudes given in the 2MASS catalog we estimated an interstellar reddening of 
E(B — V) ~ 0.04 mag. Then we estimated an average distance to the system 
as 355±7 pc. 



5. Comparison with models and other low-mass stars 

Using the radii and effective temperatures we computed the luminosities 
of the components as L x = 5.4 ± 0.4L and L 2 = 0.009 ± 0.002L for the 
primary and secondary, respectively. In Fig. 5 we compare the positions of the 
stars in the L — T P ff diagram. Th e isochrones for 5 , 6, a nd 7 Gyr obtained 
by Y 2 models f flYi et all boOjJ), f lDemarque et al.l . bflOJ)) for Z=0.004 are 
also plotted. The primary star appears to an evolved F-star with an age of 
about 6 Gyr. This comparison with the theoretical models indicates that the 
primary star should have poor me tal a bundance. 



For the first time 



Rib ad (120031 ) and iRibad (120061 ) called attention about 



the significant difference of low-mass stars' radii between measured and pre- 
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Figure 5: The primary component of T-Cygl-01385 compared with Y 2 models for Z=0.004. 
For comparison the evolutionary tracks for 1 and I.IMq and isochrones for 5, 6 and 7 Gyr 
are shown. 
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dieted by the theory. His comparison for a limited stars with masses and radii 
determined with an accuracy better than 3 percent revealed that the observed 
radii are systematically larger than those estimated by theoretical models. 
In contrary, their effective temperature s are significantly cooler than in the 



models. Later on, 



Cakirli et al. 



( 120101 ) collected the available masses, radii 
and effective temperatures of the low-mass stars obtained up to at that date 
and compared with those calculated by the stellar evolution models. They 
have also compared radii and effective temperatures of 21 low- mass stars in 
the mass-radius and mass-effective temperature diagrams. They conclude 
that stars below 0.3 M Q have radii and effective temperatures which are 
consistent with those of models. However, larger mass stars, above 0.3 M , 
begin to deviate from the theoretical predictions in the M — R diagram. The 
observed radii are significantly larger than that of models. In contrary, their 
observed effective temperatures are significantly lower than those models, i.e. 
they appear as below-shifted from the theoretical M — T e ff relation. The dis- 
crepancies in radii and effective temperatures could neither be explained by 
changing the r a tio o 



mixing-length to pressure scale height nor metallicity 
2009[ ). As an alternative explanation the magnetic activity, 
which is responsible for t he observed larger radi i but c o oler effective t emper - 



( jDemory et al. 



atures, was adopt ed by 



Ribasl (boOO l and 



Mullan k MacDonaldl (120011 ) . 



Cakirli et al. 



Torres et al. 



(120061 ) . 



(120131 ) . Due to high magnetic activity in 
fast-rotating dwarfs, because of spin-orbit synchronization, their surfaces are 
covered by dark spots. Faster rotation enhances solar-type activity which 
results in larger radius and lower effective temperature. 

Fig. 6 shows locations of the low- mass stars in the in M — R and M — T e f f 
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planes. The low-mass companion of T-Cygl-01385 appears to have a radius 
about 26 per cen t larger and effective temperature about 18 per cent cooler 



than t he models ( IBaraffe et al 



19981 ). In a very recent paper ICakirli et al. 



(120131 ) called attention to the stars having masses about 0.34 M Q . The 
empirical M — R diagram shows that the stars with masses lower than about 
0.3 M do not deviate from that predicted by the models. However, largest 
deviations both in M — R and M — T e ff occur just at the mass higher 
than 0.34 M . This mass may be taken as transition from partly convective 
atmospheres to the fully convective stars. 



6. Conclusions 

Spectroscopic observations of T-Cygl-01385 are obtained and the radial 
velocity curves of both components are revealed for the first time. Analy- 
sis of the radial velocities and the TrES at the Sloan r-passband and the 
KeplerC&m at the z-passband light curves yielded the physical parameters 
of the components. Our analysis indicates that the eclipsing binary consists 
of an F8 subgiant and an M4 dwarf. The radius of the less massive M- 
dwarf is larger about 26 percent but the effective temperature 18 
percent cooler than those estimated from the models. The radii and 
effective temperatures of low-mass stars obtained so far have been compared 
with each other and also with the models both in the M — R and M — T e ff 
panels. The less massive secondary star with a mass of O.34M ap- 
pears to be just in the boundary between partially convection and 
fully convection. 
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Figure 6: Mass-radius (left) amd mass-temperature (right) diagrams of the stars between 
0.2 to 1.0 M & in double-lined, detached ecl ipsing binaries . The c ontinuous line corresponds 



to zero-age main-sequence adopted from 



Baraffe et al 



(|1998l ) models. The radius and 



effective temperature for the less massive companion of T-Cygl-01385 are plotted as the 
open squares with standard deviations. 
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